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ABSTRACT 

Context. PAHs appear to be an ubiquitous interstellar dust component but the effects of shocks waves upon them have never been 
fully investigated. 

Aims. We study the effects of energetic (« 0.01 - 1 keV) ion (H, He and C) and electron collisions on PAHs in interstellar shock 
waves. 

Methods. We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the threshold for carbon atom loss 
from a PAH, in 50 - 200 km s~ ] shock waves in the warm intercloud medium. 

Results. Interstellar PAHs {Nc = 50) do not survive in shocks with velocities greater than 100 km s and larger PAHs (Nc = 200) 
are destroyed for shocks with velocities > 125 km s" 1 . For shocks in the « 75 - 100 km s -1 range, where destruction is not complete, 
the PAH structure is likely to be severely denatured by the loss of an important fraction (20 - 40%) of the carbon atoms. We derive 
typical PAH lifetimes of the order of a few xlO 8 yr for the Galaxy. These results are robust and independent of the uncertainties in 
some key parameters that have yet to be well-determined experimentally. 

Conclusions. The observation of PAH emission in shock regions implies that that emission either arises outside the shocked region or 
that those regions entrain denser clumps that, unless they are completely ablated and eroded in the shocked gas, allow dust and PAHs 
to survive in extreme environments. 
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1. Introduction 

Interstellar Polycyclic Aromatic Hydrocarbon molecules (PAHs) 
are an ubiquitous component of the interstellar medium. The 
mid-infrared spectrum of the general diffuse interstellar medium 
as well as energetic environments near massive stars such as H n 
regions and reflection nebulae are dominated by broad emis- 
sion features at 3.3, 6.2, 7.7, and 11.2 fim. These emission 
features are now generally attributed to infrared fluorescence 
by large PAH molecules cont aining 50-100 C-atoms, pumped 
by single FUV photons (see iTielensI 120081 for a recent re- 
view). The observed spectra also show e vidence for PAH clus- 
ters containing a few hundred C-atoms (iBregman et alJ [19891: 
Rapacioli et alj 120051 iBerne et all I2007T) as well as very small 
dust grains (~30 A; Deser t et ail 1 19901) . It seems that the in- 
terstellar grain siz e distribution extends all the way into the 
molecular domain dAllamandola et al.lfl989t iDesert et al.lll990t 
iDraine & Lil200ll) . The origin and evolution of interstellar PAHs 
are somewhat controversial. On the one hand, based upon ex- 
tensive laboratory studies of soot formation in terrestrial envi- 
ronments, detailed models have been made for the formation 
of PAHs in the ejecta of C-ri ch giants (Frenklach & Fei gelsor] 
119891: ICherchneffet"aT]ll992h - as intermediaries or as side- 
products of the soot-formation process - and studies have sug- 
geste d that such ob jects might produce enough PAHs to seed the 
ISM (lLatten fT991). On the other hand, models have been devel- 
oped where PAHs (as well as very small grains) are the byprod- 
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uct of the grinding-down process of large carbonaceous grains 
in strong s upernova shock waves which permeate the in terstellar 
medium (Borkowski& Dweklll995l: iJones etal.l[l996h . Grain- 
grain collisions shatter fast moving dust grains into small frag- 
ments and, for graphitic progenitor grains, these fragments might 
be more properly considered PAH molecules. The destruction 
of interstellar PAHs is equally clouded. Laboratory studies have 
shown that small (less than 16 C-atoms), (catacon densed) P AHs 
are ra pidly photodissociated by ~ 10 eV photons ( Jochims et al. 
1994). However, this process is strongly size-dependent as larger 
PAHs have many more modes over which the internal energy can 
be divided and PAHs as large as 50 C-ato ms might actually be 
stable against photodisso ciation in the ISM (Le Page et al. 2001; 
lAllamandola et al] Q989). While strong shock waves have been 
considered as formation sites for interstellar PAHs, the destruc- 
tion of these PAHs in the hot postshock gas has not been eval- 
uated. Yet, high energy (~ 1 keV) collisions of PAHs with ions 
and electrons are highly destructive. 

The observational evidence for PAHs in shocked regions is 
quite ambiguous. The majority of supernova remna nts does not 
show PAH feature s (e.g. Cas A. ISmifh et alJl2008h . but obser- 
vations of N132D dTappe et al.ll2006l) su ggest the possib i lity of 
PAH survival in shocks. Recent work by Andersen et al. (2007) 
investigates the presence of PAHs in a subset of galactic super- 
nova remnants detected in the GLIMPSE survey. Unfortunately 
the interpretation of such observations is not straightforward, be- 
cause of the difficulty in disentangling the PAH features intrin- 
sic to the shocked region with those arising from the surround- 
ing material. Another interesting case is the starburst galaxy 
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M82, which shows above and below the galactic plane a huge 
bipola r outflow of shock- heated gas interwoven with PAH emis- 
siorfl dArmus et al.1 l2007). PAHs have also been observed at high 
galact ic latitudes in the edge-on galaxies NGC 5907 and NGC 
5529 dlrwin & Madden! 120061: llrwin et al.ll2007l) . Shock driven 
winds and super novae can create a so-called "galactic fountain" 
(Bre gmanlll980t) transporting material into the halo and these 
detections of PAHs suggests the possibility of survival or for- 
mati on of the molecul e s under thos e conditions. On the other 
hand lO'Halloran et all (120061 120081) have found a strong anti- 
correlation between the ratio [FeII]/[NeII] and PAH strenght in 
a sample of low-metallicity starburst galaxies. Since [Fell] has 
been linked primarly to supernova shocks, the authors attributed 
the observed trend to an enhanced supernova activity which led 

to PAH destruction. 

In our previous study dJones et al.| [T996). we considered the 
dynamics and processing of small carbon grains with Nc > 100. 
The processing of these grains by sputtering (inertial and ther- 
mal) in ion-grain collisions and by vaporisation and shattering in 
grain-grain collisions was taken into account for all the consid- 
ered grain sizes. In that work, the smallest fragments (af < 5 A) 
were collected in the smallest size bin and not processed. In this 
work we now consider what happens to these smallest carbon 
grain fragments that we will here consider as PAHs. In this paper, 
we will consider relatively low velocity (< 200 km s~') shocks 
where the gas cools rapidly behind the shock front but, because 
of their inertia, PAHs (and grains) will have high velocity col- 
lisions even at large postshock column densities. Collisions be- 
tween PAHs and the gas ions occur then at the PAH velocity 
which will slowly decrease behind the shock front due to the 
gas drag. This relative velocity is thus independent of the ion 
mass and, for dust grains, destruction is commonly called iner- 
tial sputtering. Destruction of PAHs in high velocity (> 200 km 
s _1 ) shocks - which cool slowly through adiabatic expansion - 
is dominated by thermal sputt ering and these shock s are con- 
sidered in a subsequent paper dMicelotta et al.l 120 09, hereafter 
MJT). 

This paper is organized as follows: Sect. 2 describes the the- 
ory of ion interaction with solids, Sect. 3 illustrates the appli- 
cation of this theory to PAH processing by shocks and Sect. 4 
presents our results on PAH destruction. The PAH lifetime in 
shocks and the astrophysical implications are discussed in Sect. 
5 and our conclusions summarized in Sect. 6. 



2. Ion interaction with solids 

2. 1 . Nuclear interaction 

The approach used in our earlier work is not valid for planar 
PAH molecules with of the order of tens of carbon atoms. Here, 
we assume that co llisions are binary in nature, as is assume d 
in work on solids dLindhard et al.lll963l Il968t ISigmundlll98lh . 
If the energy transfer is above the appropriate threshold value, 
we assume that the carbon target is ejected from the molecule. 
For energy transfer below that threshold, the energy will become 
thermal energy and be radiated away. 

In this description the "bulk" nature of the target enters only 
after the first interaction, when the projectile propagates into the 
material. We therefore consider only the first interaction, which 
is described in the binary collision approximation in a way that 
then conveniently allows us to take into account the "molecular" 
nature of the target. 

1 http://chandra.harvard.edu/photo/2006/m82/ 



In addition to the energy directly transferred to the target nu- 
cleus through elastic scattering (nuclear stopping or elastic en- 
ergy loss), the energy loss to the atomic electrons (electronic 
stopping or inelastic energy loss) should also be considered 
(Lindh ard etal.1 [19631 [ 1968). In a solid the energy transferred 
via electronic excitation is distributed around the impact region. 
For a PAH, which has a finite size, the energy will be spread 
out over the entire molecule. This energy will either be radiated 
away or a fragment can be ejected. 

Nuclear and electronic stoppin g are simultaneous p rocesses 
which can be treated separately (Lind hard et al.|[l963l) . Fig. QT| 
illustrates these effects and shows the PAH evolution following 
the loss of carbon atoms, /Vc(lost), for the two limiting cases: 
1) where there is an instantaneous and random removal of the 
lost carbon atoms and 2) where the carbon atoms are removed 
only from the periphery in order to preserve aromatic domain as 
much as possible. The reality of PAH erosion in shocks prob- 
ably lies somewhere between these two extremes and will in- 
volve isomerisation and the formation of five-fold carbon rings 
that distort the structure from a perfectly two-dimensional form. 
This then begs the question as to the exact form and structure of 
small carbon species once growth resumes by atom insertion and 
addition. The full treatment of the nuclear stopping is given here, 
for the electronic stopping only the results of the calculations are 
shown, for the complete description of the phenomenon we refer 
the reader to paper MJT. 

The treatment of PAH processing by shocks should also in- 
clude the effects of fast electrons present in the gas. Because 
of their low mass, electrons can reach high velocities and 
hence high collision rates even at relatively low temperatures 
(T ~10 5 K), leading to potentially destructive collisions. Again 
for a detailed description of the electron-PAH interaction see pa- 
per MJT. 

The theory of ion penetration into solids described here con- 
siders collisions where the transferred energy T goes from to 
the maximum transferable energy. For this study, we are inter- 
ested in only those collisions that are able to remove carbon 
atoms from the PAH, i.e. for which the energy transferred is 
greater than the minimum energy Tq required for C ejection. In 
Sect. 2.2 we present the modifications we introduce into the the- 
ory in order to treat the case of collisions above this threshold. 

To describe the binary collision betwe en a moving at om 
(or ion) and a stationary target atom (e.g. lSigmund] ll98ll) . a 
pure classical two-particle model using the Coulomb repulsion 
between the nuclei (Rutherford scattering) is adequate only at 
high energies, i.e. when s » 1, where s is the dimensionless 
Lindhard' s reduced energy 



Mi + M 2 Zi Z 2 e 2 

where M\ and Z\ are the mass and atomic number of incident 
particle respectively, M 2 and Z2 the mass and atomic number of 
target particle, E is the kinetic energy of incident particle and e 
is the electron charge, with e 1 - 14.39 eVA. The quantity a is the 
screening length, a parameter that defines the radial spread of the 
electronic charge about the nucleus. For the screening length we 
adopt the Univ ersal Ziegler - Bier sack - Littmark (ZBL) screen- 
ing length a\j (IZiegler et al.H l 985) 

flu = 0.885 fl (Z? 23 + Z" 23 )- 1 (2) 

where ao = 0.529 A is the Bohr radius. The condition s » 1 
implies that the energies are large enough that the nuclei ap- 
proach closer to each other than the screening length a. At lower 
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LABORATORY COORDINATES 

Fig. 1. Scattering geometry for an elastic collision of particle 1 
(mass Mi , initial velocity vq, impact parameter p), on particle 2 
(mass M 2 , initial velocity zero). After the impact, the projectile 
particle 1 is deflected by the angle # and continues its trajectory 
with velocity vi. The target particle 2 recoils at an angle (f> with 
velocity \ 2 , 

Table 1. Kinetic energy E and reduced energy e for H, He and 
C impacting on a carbon atom. 



Projectile 


£l(v pl ) 




E a (v p2 ) 


e*(v p2 ) 


H 


7.30 


0.031 


117.4 


0.50 


He 


29.4 


0.048 


469.7 


0.76 


C 


88.1 


0.028 


1409. 


0.45 



" Kinetic energy in eV. The projectile velocities v p i = 37.5 km s and 
v p2 = 150 km s -1 are defined by the shock velocity v$ via the equation 



Vs, with vs = 50 and 200 km s 1 respectively. 
Dimensionless, calculated from Eq.fT]and Eq.|2] 



energies, s < 1, it is essential to consider the screening of the 
Coulomb interaction. In this case the Rutherford approximation 
is not adequate and the scattering problem must be treated using 
a different approach. 

To choose the appropriate formalism to describe our interac- 
tion, we need to calculate the reduced energy for our projectiles. 
For our study of the behaviour of PAHs in shocks, we consider 
the binary collision between H, He and C ions (projectiles) and 
a carbon atom (target) in the PAH molecule. The velocity v p of 
the projectile is determined by the shock velocity vs through the 
relation v p = | vs- We consider here shock velocities between 50 
and 200 km s _1 . The corresponding projectile kinetic energies E 
and reduced energies e are reported in Table Q] for the two limit- 
ing cases v p i = |(50) = 37.5 km s~' and v P 2 = f (200) = 150 km 
0-1 



The calculation clearly shows that for the shocks we are con- 
sidering e < 1, implying that our problem cannot be treated in 
terms of Rut herford scattering but requires a different formalism, 
described by Sigmund (1981) and summarized below. The scat- 
tering geometry for an elastic collision of the projectile particle 
1 on target particle 2 is illustrated in Fig.Q] Particle 1 has mass 
Mi, initial velocity vo and impact parameter p, where the impact 
parameter is the distance of closest approach of the centers of the 



Table 2. Threshold energy To and critical kinetic energy Eq d for 
H, He and C ions impacting on a carbon atom. 



£on(H) £bn(He) £ 0n (C) 



4.5 


15.8 


6.0 


4.5 


7.5 


26.4 


10. 


7.5 


10. 


35.2 


13. 


10. 


12. 


42.3 


16. 


12. 


15. 


52.8 


20. 


15. 



Note: Tn and En„ in eV. 



two atoms/ions that would result if the projectile trajectory was 
undeflected. Particle 2 has mass M 2 and is initially at rest. After 
the impact, the projectile is deflected by the angle § and contin- 
ues its trajectory with velocity vi. A certain amount of energy 
T is transferred to the target particle which recoils at an angle <p 
with velocity v 2 . The maximum transferable energy corresponds 
to a head-on collision (impact parameter p — 0) and is given by 



y E 



\M X M 2 
(Mi + M 2 ) 2 



(3) 



An important quantity to consider is the nuclear stopping 
cross section S n (E), which is related to the average energy loss 
per unit path length of a p article travelling throu gh a material of 
atomic number density N (Lindh ard et al Jl 19631) 



dE 
dR 



J' 



NS n (E) = N AaiE, T) ■ T 



(4) 



where cr(E, T) is the energy transfer cross section (see Appendix 
A for details). S n {E) has the dimensions of (energy x area x 
atom -1 ) and in fact represents the average energy transferred per 
atom in elastic collisions when summed over all impact param- 
eters. 

The nuclear stopping cross section can be expressed in terms 
of the Lindhard's reduced energy s and th e dimensionless re- 
duced nuclear stopping cross section s a (s) dLindhard et ali i 968. 
see Eq. IA.4I and \A.5i . For this latter we adopt the Universal 
reduced Ziegler - Biersack - Littm ark (ZBL) nuclear stopping 
cross section dZiegler et al.ll 19851) . which is an analytical ap- 
proximation to a numerical solution that reproduces well the ex- 
perimental data. The ZBL reduced nuclear stopping cross sec- 
tion has the form 



0.5 ln(l + 1.1383 s) 



e + 0.01321s' 

In s 
~2e~ 



0.21226 . 



0.19593 s 1 



.0.5 



e< 30 



e> 30 



(5) 



and the nuclear stopping cross section S n (E) can be written as 



S n {E) — 4 n flu Z\ Z2 e 



Mi +M 2 



(6) 



with the screening length flu from Eq. [2] 
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2.2. Nuclear interaction above threshold 

For this study we are interested in destructive collisions, i.e., col- 
lisions for which the average transferred energy T exceeds the 
minimum energy To required to remove a carbon atom from the 
PAH. The theory discussed in Sect. 2.1 does not treat this situa- 
tion and considers the specific case where To = (no threshold). 
To include the treatment of collisions above threshold (To > 0) 
we developped the appropriate expressions for the relevant quan- 
tities described in the previous sections. 

The definition of the nuclear stopping cross section S n (E) 
can be written in a more general way as 



Sn(E) 



-r 



do-(E, T) ■ T 



(7) 



where Tq > 0. The total energy transfer cross section per carbon 
atom cr(E) is defined by 



o-(E) 



-P 

J T 



do-(E, T) 



(8) 



In this case the threshold To must be strictly positive, otherwise 
cr would diverge (this can be verified by substituting the expres- 
sion for dcr from Eq. I A. 1 1 and evaluating the integral). Finally, 
the average energy transferred in a binary collision is given by 
the ratio between S „ and cr 



(T(E)} 



o-(E) 



(9) 



The condition yE = T m > To — y Eq„ then imposes, for the 
kinetic energy of the incoming ion, the condition that E > Eo n = 
Toly. In Table|2]we report the critical energies £o n for H, He and 
C ions corresponding to different values of the threshold energy 
7"o. 

Using dcr from Eq. I A. 1 1 and evaluating the above integrals 
we obtain 



S a (E) = 



C F~ m 
1 — m 



ryl-m rp I— mi 
l 1 m 1 Q J 



1 —m i 



o-(E) 



C E 



(T(E)) = 



m J 



j^l—m 
1 



r r— m 

'o 



(10) 



(ID 



(12) 



To calculate the quantity m = m(E) we use the following 
expression from lZiegler etafl (fl985) 



m(E) — 1 - exp 



- exp 



2 at (o.l In ( 

!=0 V ^ 



s(E) 

El 



(13) 



with s } = 10 } and a-, = -2.432, -0.1509, 2.648, -2.742, 1.215, 
-0.1665. 

Combining Eq. IA.3I and Eq. IA.41 after some algebraic ma- 
nipulation, we can rewrite the above expressions for S n , cr and 
(T) in the more convenient form shown below. The full calcu- 
lation is reported in Appendix A. As explained in Sect. 2.1, we 
adopt for the reduced stopping cross section the ZBL function 
■s„ (e) (Eq.|5]l with the appropriate screening length a\j. 



cr(E) = AnaZxZze^ 



Mi 



Mi + M 2 



1 — m 1 



y E 



(tr- 



(T(E)} 



pi— m jj\—tn 

11 ^On 



1-m' E^-E- 



(15) 



(16) 



Note that the term outside of the square brackets in Eq.[T4land[T5l 
is the stopping cross section S n (E) when To = (no threshold). 

The nuclear stopping cross section S n (E) (eV A 2 atom -1 ), 
the total energy transfer cross section cr(E) (A 2 atom -1 ) and the 
average energy transferred (T{E)) (eV), for H, He and C ions im- 
pacting on a carbon atom, calculated from the above expressions 
assuming a threshold Tq = 7.5 eV, are shown in Fig. [2] 

The sharp cut on the left-hand side of the curves arises from 
the fact that we are treating collisions above threshold, and these 
quantities are defined only for energies of the incident ion greater 
than the critical value £on- It can be seen that all quantities in- 
crease in absolute value with increasing atomic number and mass 
of the projectile (Z\ and Mi). The two vertical lines indicate the 
minimum and maximum kinetic energy of the projectile consid- 
ered in our study, corresponding to the PAH velocity in the 50 
and 200 km s _1 shocks respectively. The values are those cal- 
culated in Sect. 2.1 and reported in Table [1] The figure cleary 
shows that for hydrogen the critical value Eqo is greater than 
the lower limit of energy range. This implies that in the lower 
velocity shocks hydrogen is not energetic enough to cause car- 
bon ejection. The curves for S n presents a characteristic convex 
shape with a maximum, illustrating that nuclear energy transfer 
is important only for projectiles with energy falling in a specific 
range. In particular, the nuclear stopping becomes zero at high 
energies, with a limiting value depending on projectile and tar- 
get: in our case, going from H to C impacting on carbon, the 
curves extend further to the right, in the direction of higher ener- 
gies. In the high energy regime, the energy transfer is dominated 
by electronic stopping (MJT). 

For a given incident ion energy, the difference between the 
values of S n in the threshold and no-threshold cases results from 
the definition of the nuclear stopping and from the properties of 
dcr. The differential cross section (cf. Eq. lA.U strongly prefers 
collisions with low energy transfers (T <K T m ) and, moreover, 
decreases in absolute magnitute with increasing E. For each E, 
S n is defined as the integral over the transferred energy T, of 
the product between T and the corresponding cross section dcr. 
Choosing To > means excluding from the integral all energy 
transfers T < To, for which the cross section has the highest 
values. The remaining terms have higher values of T but lower 
values of dcr, then the integral gives a result smaller than the no- 
threshold case, which includes all small energy transfers with 
their higher cross sections. 

The total cross section cr(E) clearly shows that the projec- 
tiles can efficiently transfer energy to the target atom only when 
their kinetic energy lies in the appropriate window. In particular, 
it can be seen that the average energy transferred (T) increases 
with E, nevertheless at high energies cr is close to zero (and the 
collision rate will be small). For a fixed target atom (in our case, 
carbon), the width of the cr curve, and consequently the width 
of the energy window, increases with Z and M of the projectile. 
Heavier ions transfer more energy and in a more efficient way. 



S n (E) = AnaZiZie 1 



Mi 

Mi + M 2 



(17 \ 1-m 
-f) 



(14) 
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Energy Incident Ion (eV) Energy Incident Ion (eV) 

Fig. 2. The nuclear stopping cross section S n (E), the total cross section cr(E) and the average energy transferred (T(E)) calculated 
for H, He and C ions impacting on a carbon atom. The curves are calculated for the threshold energy 7q = 7.5 eV. The nuclear 
stopping cross section 5° corresponding to To = (no threshold) is shown for comparison. The two vertical lines indicate the 
limiting energies for the incident ion. These are defined as the kinetic energies of the projectile when its velocity v p equals | (vsi,s2), 
where vsi = 50 km s" 1 and vs2 = 200 km s _1 are the lowest and highest shock velocities considered in this study. 
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2.2.1. The threshold energy T 

The threshold energy, 7b, is the minimum energy that must be 
transferred via nuclear excitation to a carbon atom, in order to 
eject that same atom from the PAH molecule. The choice for 
To for a PAH is unfortunately not well-constrained. There are 
no experimental determinations, and the theoretical evaluation is 
uncertain. The analog of To in a solid is the displacement energy 
Td, defined as the minimum energy that one atom in the lattice 
must receive in order to be moved more than one atomic spac- 
ing away from its initial position, to avoid the immediate hop 
back into the original site. For graphite, the data on the thresh- 
old energy for atomic displacement differ significant ly, varying 
from -30 eV (Montet 1967t iMontet & Mverslfl97ll) to 12 eV 
dNakai et~aT1 1 1 99 lh largely depending on direction (eg., within 
or perpendicular to the basal plane). For a PAH, the lower value 
(corresponding to the perpendicular direction) seems then more 
appropriate. For amorphous carbon, Cosslett (1978) h as found 
a low value of 5 eV. Electron microscopy studies by iBanharj 
(1 19971) on graphitic nanostructures irradiated with electrons of 
different energies, indicate that a value of 7^ ~ 15-20 eV seems 
appropriate for the perpendicular direction. The in-plane value, 
however, could be much higher, presumably above 30 eV. 

Instead of graphite, fullerenes and carbon nanotubes may be 
a better analog for PAH molecul es. For fullerene, T& h as been 
found between 7.6 and 15.7 eV dFiiller & Banhartlll996l) . Single 
walled nanotubes consist of a cylindrically curved graphene 
layer. Unfortunately, also in this case the threshold for atomic 
displacement is not precisely determined. However it is expected 
to be lower than in a multi-la yered tube, fo r which a value of 
I'd ~ 15-20 eV has been found dBanhartl 19971) close to the value 
of graphite. We note that 4.5 and 7.5 eV are close to the energy 
of the single and double C-bond. 

Because we cannot provide a well-defined Tq, we decided to 
explore a range of values, to study the impact of the threshold 
energy on the PAH processing. For our standard case, we adopt 
7.5 eV that we consider a reasonable value consistent with all 
the experimental data. However, we have varied To from 4.5 to 
15 eV. 

Fig. [3] shows the comparison between S n , cr and {T) calcu- 
lated for He on C assuming To = 4.5, 7.5 and 15 eV. Coherently 
with their definition, S „ and <x increase with decreasing thresh- 
old, because more collisions are effective and the cross section 
increases with decreasing energy. Of course, the average energy 
transferred will decrease when the threshold energy is decreased. 

In Sect. 4. 1 we discuss the effect of the choice of different 
values for Tq on the PAH survival in shocks. 



3. PAHs in shocks 

When grains and PAHs enter a shock they become charged and 
then gyrate around the compressed magnetic field lines. This 
leads to relative gas-particle velocities and hence to collisions 
with the gas (and other grains/PAHs). Collisions with the gas 
result in drag forces and therefore a decrease in the relative gas- 
particle velocity. However, these same collisions with the gas 
can also lead to the removal of atoms from the particle if the 
relative velocites are larger than the given threshold for an ero- 
sional process. The removal of carbon atoms from the PAH due 
to ion collisions, where the impact velocity is determined by the 
relative motion between the two partners, is the analog of the in- 
ertial sputtering of dust particles due to ion-grain collisions. In 
the following we will then refer to it using the term inertial, and 
the same will apply for all the related quantities. 
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Energy Incident Ion (eV) 

Fig. 3. The nuclear stopping cross section S n {E), the cross sec- 
tion cr(E) and the average energy transferred (T(E)) calculated 
for He ions impacting on a carbon atom, calculated for three 
values of the threshold energy Tq'. 4.5, 7.5 and 15 eV. The two 
vertical lines delimit the energy range of interest (cf. Fig. 



In determining the processing of PAHs in shock waves, 
as with all grain processing, it is the relative gas-grain veloc- 
ity profile through the shock that determines the level of pro- 
cessing. In calculating the relative ion-PAH velocity through 
the shock we use the sa me approach as in our previous work 
dJones et aL[ 1994 1 19961). whi ch is based on the methods de- 
scribed in lMcKee et al.l dl987l) . The PAH velocity is calculated 
using a 3D particle of the same mass as the 50 carbon atom PAH 
under consideration. The PAH velocity depends then on the PAH 
mass and average geometric cross section. For a PAH with Nc 



carbon atoms, these are given by Ncmc and 0.57rap AH with apah 
given by 0.9 y/Nc A, appropriate for a compact PAH dOmontl 
1986) and the factor 1 /2 in the cross section takes the averaging 
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Fig. 4. The structure of the 125 km s _1 shock: temperature T\ = 
77 10 4 K, hydrogen density «h and electron relative abundance 
X e - All quantities are plotted as a function of the shocked column 
density Ah = «o vs t, where the preshock density «o = 0.25 cm" 3 . 
To convert column density to tim e, use the follo wing relation: 
log f(yr) = log A H (cnr 2 ) - 13.9 dJones et al.lll996h . 
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Fig. 5. The velocity profile of a 50 C-atom PAH in a shock with 
velocity vs = 125 km s _1 . The PAH velocity vpah, plotted as a 
function of the shocked column density Ah, represents the rel- 
ative velocity between the molecule and the ions present in the 
shock. Overlaid is the PAH charge Zpah along the sho ck, calcu- 
lated using the theory described in McK ee et al.l (Il987h . 



over impact angle into account. The PAH and grain cross sec- 
tions are very close (to within 1 1 % for Ac = 50), thus we are 
justified in using the same numerical approach even though we 
are using a 3D grain to calculate the velocity profile of a 2D PAH 
through the shock. 

The PAHs are injected into the shock with 3/4 of the shock 
speed, as are all grains, and their trajectories are then calculated 
self-consistently with their coupling to the gas, until the rela- 
tive gas-PAH veocity becomes zero. The velocity calculation in- 
cludes the effects of the direct drag with the gas due to atom 
and ion collisions and the drag due to the ion-charged PAH in- 
teraction in the post-shock plasma. We find that for some shock 
velocities, in our case for vs = 75 and lOOkms" 1 , the PAHs 
(and grains) experience betatron acceleration in the post-shock 
gas. All the relevant expressions and assumptions for the cal- 
culation of the grain velocit y, betatron accelera tion and grain 
charge are fully described in Mc Kee et alJ dl987l) . Thus, in cal- 



Fig. 6. The number of collisions A t = /?„j x t of a 50 C-atom 
PAH with H, He and C ions in a 125 km s _1 shock, as a function 
of the shocked column density, Ah- The cross section has been 
evaluated for a threshold energy, To, of 7.5 eV. The carbon curve 
has been multiplyied by a factor 100 for comparison. 



culating the post-shock PAH velocity profiles, we follow exactly 
the same methods as used in our previous work. The structure of 
the 125 km s" 1 shock is shown in Fig.|4]as a function of the col- 
umn density Ah- Fig|5]shows the velocity profile for a 50 carbon 
atom PAH in the same shock, together with the effective charge 
of the molecule, used to calculate the velocity profile itself. The 
50 C-atoms PAH is positively charged (charge between +2 and 
+3) during the whole slowing process, and approaches neutrality 
at the end of the shock. 



3.1. Ion collisions: nuclear interaction 

Knowing the velocity profile of the PAH, we can then calculate 
the inertial collision rate PAH-ions R n \ (s ) through the shock. 
This is given by the following equation 



R n ,i(Nn) = 0.5 Xi »h vpah o" A c F c 



(17) 



where «h(A^h) is the hydrogen particle density along the shock, 
Vpah(A^h) the PAH-ion relative velocity along the shock and 
Xi is the relative abundance of the projectile ion with respect 
to hydrogen. We adopt the gas phase abundances xu ■ XKe '■ 
Xc — 1 : 10" 1 : 10~ 4 , where the carbon abun dance is bet ween 
the values (0 5 - l)x l0~ 4 and 1.4xl0~ 4 from [Sofia] d2009b and 
ICardelli etal] (1 19961) respectively. 

The term ct(Nh) is the cross section averaged over those col- 
lisions that transfer an energy larger than the threshold energy To 
per C-atom and this cross section should therefore be multiplied 
by the number of carbon atoms in the PAH, Ac- The factor 0.5 
takes the angle averaged orientation into account (see appendix 
C). 

Because both collision partners are charged, the effect of 
the Coulombian potential must be included as well. Depending 
on whether the interaction is attractive or repulsive, the en- 
ergy transfer cross section will be increased or reduced by the 
coulombian factor Fc given by 



F c = 1 



2 Zi on ZpAH e 



4 n £q apAH Mi m H v^ AH 



(18) 



where Z; on = +1 and M\ are the charge and the atomic mass 
(in amu) of the incident ion, Zpah is the charge of the PAHs 
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along the shock, «pah is the PAH radius and Vpah the PAH-ion 
relative velocity. The constant e is the electron charge, sq is the 
permittivity of the free space and toh the mass of the proton. 
The total number of destructive collisions is then given by the 
integral of the collision rate (Eq. Wl\ behind the shock, 

N t = J RniiNn) d? (19) 

where it should be understood that the postshock column density 
A^h and the time, t, are related through A^h = no^s?- With the 
proper cross section, this number A^ t , is then equal to the number 
of carbon atoms lost by a PAH in collisions with H, He, or C. 

Figure|6]illustrates the destructive collisions for a 50 C-atom 
PAH behind a 125 km s _1 shock assuming To = 7.5 eV. These 
results are plotted in such a way that equal areas under the curve 
indicate equal contributions to the total number of destructive 
collisions. R n j drops precipitously because of the drop in rela- 
tive PAH-gas velocity. Because heavier projectiles are more en- 
ergetic in inertial collisions, this drop off shifts to higher col- 
umn densities for heavier species. The results show that He is 
much more effective in destroying PAHs than H because of the 
increased energy transferred for heavier collision partners (cf. 
Figure [2]). The low abundance of C depresses its importance in 
inertial sputtering. 

The number of carbon atoms in a PAH is now given by 

N c (0 = N c (0) exp [-NJNc] (20) 

and the fraction of carbon atoms ejected from this PAH is 

F L = (1 - expHVt/Afc]) (21) 

where N t is now evaluated throughout the shock. 

In the shocked gas, the velocity of the ions is not only de- 
termined by the relative motion with respect to the PAH (inertial 
case), but also by the temperature of the shocked gas. In princi- 
ple, the inertial and thermal velocity should be added vectorally 
and averaged over the angle between the inertial motion and the 
(random) thermal motion as well as over the thermal velocity 
distribution. However, that becomes a quite cumbersome calcu- 
lation and, hence, we will follow c alculations for spu ttering of 
dust grains in interstellar shocks Ccf. lJones et alJ fl994) and eval- 
uate these two processes (inertial and thermal sputtering) inde- 
pendently. Studies have shown that this reproduc es more exten- 
sive calculations satisfactorily dGuillet et al.ll2.Q07b . The thermal 
destruction rate is given by 

R n j(N H ) = Nc0.5xinn F c (v) vcr(v)/(v, T) dv (22) 
Jvo 

with / (v, T) the Maxwellian velocity distribution. The tempera- 
ture has to be evaluated along the shock profile (cf. Figure[4]i and 
care should be taken to only include velocities corresponding to 
energies larger than the threshold energy, To (e.g., with E > Eo n ; 
cf., Table[2]). The fraction of C-atoms ejected by this process can 
be evaluated analogously to Eq. [21] 

3.2. Ion collisions: electronic interaction 

As reported in the introduction of the paper, the collision be- 
tween PAH and ions triggers two simultaneous process, which 
can be treated separately: the nuclear stopping (elastic energy 
loss) and the electronic stopping (inelastic energy loss). The first 
has been extensively discussed in the previous sections, while 
for the full treatment of the electronic interaction we refer the 



reader to MJT For the sake of clarity, we report here the essen- 
tial concepts and the principal equations which will be used in 
the following. 

The energy transferred to the electrons is spread out over 
the entire molecule, leaving the PAH in an excited state. De- 
excitation occurs through two pricipal decay channels: emission 
of infrared photons and dissociation and loss of a C2 fragment. 
This latter is the process we are interested in, because it leads to 
the PAH fragmentation. The dissociation probability p (see Sect. 
4. 1 in MJT) depends on the binding energy of the fragment Eq, 
on the PAH size, Nc, and on the energy transferred, which in 
turns depends on the initial energy (velocity) of the projectile. 

For a fixed value of the transferred energy, the dissociation 
probability decreases for increasing Eq and Nc because either 
more energy is required in the bond that has to be broken or 
because the energy is spread over more vibrational modes and 
hence the internal excitation temperature is lower. On the other 
hand, the more energy that is deposited in the PAH, the higher 
is the dissociation probability. The energy transferred via elec- 
tronic excitation increases with the energy of the projectile up to 
a maximum value, corresponding to an incident energy of 100 
keV for H (and higher for more massive particles), and decreases 
for higher energies. The deposited energy also increases with the 
path-length through the molecule and will be higher for larger 
PAHs impacted at grazing collision angles. For the shocks con- 
sidered in this study, the energy transferred increases with in- 
cident energy (velocity) and hence the dissociation probability 
increases as well. 

As for the nuclear stopping, also for the electronic interac- 
tion we have to consider the effect of both inertial and thermal 
velocities. The inertial collision rate is given by 

R s ,i(Nn)=vvmXinnF c o- g (#) />(v P ah, #) sintfdtf (23) 

Ji?=0 

where § is the angle between the axis normal to the PAH plane 
and the direction of the incoming ion. The term cr g is the geo- 
metrical cross section seen by an incident particle with direction 
defined by §. The PAH is modelled as a thick disk with radius 
«pah and thickness d, then the cross section is given by 

cr g = n ap AH cos § + 2 ami! d sin # (24) 

which reduces to <x g = 7rapA H for § = (face-on impact) and to 
<x g = 2apAH^ for § - n/2 (edge-on impact). The term p(vp^H, &) 
represents the total probability for dissociation upon collision 
via electronic excitation, for a particle with relative velocity vpah 
and incoming direction § (see Sect. 4.1 in MJT). 
For the thermal collision rate we have 

R e j(N H ) = R eJ (v)f(v,T)dv (25) 
Jvo 

where the temperature T = T(Nu) is evaluated along the shock. 
The lower integration limit vo is the ion velocity corresponding 
to Eq. The number of carbon atoms lost can be evaluated anal- 
ogously to the nuclear interaction but care should be taken to 
include the loss of 2 C-atoms per collision. 

There is a clear distinction between the nuclear and elec- 
tronic interactions. In nuclear interactions, a C-atom is ejected 
because a direct collision with the impacting ion transfers 
enough energy and momentum to kick out the impactee instan- 
taneously. In electronic interaction, the impacting ion excites the 
electrons of the PAH. Internal conversion transfers this energy 
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to the vibrational motions of the atoms of the PAH. Rapid in- 
tramolecular vibrational relaxation leads then to a thermaliza- 
tion of this excess energy among all the vibrational modes and 
this can ultimately lead to dissociation (or relaxation through IR 
emission). The threshold energy in the nuclear process, To, dif- 
fers therefore from the electronic dissociation energy, Eq. The 
latter really is a parameter describing the dissociation rate of a 
highly excited PAH molecules using an Arrhenius law and this 
does not necessarily reflect th e actual bindi ng energy of the frag- 
ment to the PAH species (cf. lTielensl[2005h . Following MJT, we 
will adopt the canonical value of 4.6 eV for Eq. However, this 
energy is very uncertain and we will evaluate the effects of re- 
ducing and increasing this parameter to a value of 3.65 and 5.6 
eV respectively (MJT). 

3.3. Electron collisions 

For the full treatment of the PAH collisions with electrons, we 
refer again to the paper MJT, providing here a short summary of 
the basic concepts and equations. 

Because of their small mass, the thermal velocity of the elec- 
trons always exceeds the inertial velocity of the PAH. Hence, 
only the thermal destruction needs to be evaluated. We follow 
the same formalism used for the electronic interaction in ion- 
PAH collisions. The energy dumped into the molecule during 
collisions with electrons is spread over and determines (with Eq 
and Nc) the value of the dissociation probability. The electron 
energy loss rises sharply with the electron energy, reaching its 
maximum for incident energy around 100 eV. This energy range 
falls exactly in the interval relevant for our shocks, implying that 
the electrons optimally transfer their energy. 

The thermal electron collision rate can be written as 

flelec,T(jV H ) = I S(v)/(V,T)dV (26) 

J v '0,clec 

S(v) = vx e n H F Ce o- g (&)p(v,&)sin§d& (27) 
J#=o 

where vo e i ec is the electron velocity corresponding to Eq and Xe 
is the electron relative abundance along the shock. The electron 
coulombian factor Fc e is always equal to 1 (within less than 1%) 
because electrons have low mass and high velocities with respect 
to ions. The temperature T evaluated along the shock is the same 
as for the ions, but electrons will reach much larger velocities. 
From Eq.l26landl27lwe expect then to find a significantly higher 
collision rate with respect to the ion case. The fraction of C- 
atoms lost by electron collisions can be evaluated analogously 
to that for ions (cf. Eq. l2Tb. 

4. Results 

Fig.|7]and[8]show the fraction of carbon atoms ejected from a 50 
and 200 C-atoms PAH due to collisions with electrons and H, He 
and C, assuming the nuclear threshold energy To = 7.5 eV and 
the fragment binding energy Eq = 4.58 eV. The results concern- 
ing nuclear, electronic and electron interaction are discussed in 
the following sections. 

4. 1 . PAH destruction via nuclear interactions 

For the inertial nuclear interactions, the fraction of ejected car- 
bon atoms Fl depends on both cr and x- Hydrogen has the high- 
est abundance (xn - 1) but the lowest absolute value for the cross 



section (see Fig. |2j. In addition, cr is significantly different from 
zero only for the highest shock velocities. This results in contri- 
bution to atom ejection which is only relevant for vs above 150 
km s _1 . Helium is ten times less abundant than hydrogen (xm = 
0. 1), but this is compensated for by a higher cross section for all 
shock velocities. In particular, the C-atom ejection curve shows a 
peak between 50 and 125 km/s due to betatron acceleration: be- 
cause of the higher velocity, the collision rate increases (cf. Eq. 
[PTT i and then the PAHs experience more destructive collisions. 
After the peak, as expected the curve increases with the shock 
velocity. In the case of carbon, the increased cross section is not 
sufficient to compensate for the low abundance (xc = 10~ 4 ), re- 
sulting in a totally negligible contribution to PAH destruction. 
For all shock velocities, the fraction of C-atoms removed be- 
cause of inertial nuclear interaction does not exceed the value of 
20%. 

Concerning the thermal nuclear interaction, carbon does not 
contribute to PAH destruction because of its very low abundance 
compared to H and He, as for the inertial case. For hydrogen and 
helium, as expected for low velocity shocks, the temperature is 
generally not sufficiently high to provide the ions with the en- 
ergy required to remove C-atoms. Nevertheless, the ions in the 
high velocity tail of the Maxwellian distribution can be energetic 
enough to cause C-atom ejection, as can be seen for He at 100 
km s _1 . This is less evident for hydrogen. In this case the criti- 
cal energy £on is higher than for helium and carbon. The corre- 
sponding critical velocity vo will be higher as well. For the lower 
velocity shocks, the peak of the hydrogen maxwellian function 
f(v, T) is well below Vo, as a consequence the integrand of Eq. 
l22lis close to zero over the integration range, and the same will 
be true for the collision rate. At the highest shock velocities the 
curves show a similar trend, with a steep rise beyond 125 km 
s _1 leading to complete PAH destruction, i.e. removal of ALL 
carbon atoms, for shock velocities above 150 km s _1 . At around 
135 km s the hydrogen contribution becomes larger than that 
for helium. At these high velocities the He and H cross sections 
reach approximately their maximum values (cf. Fig. |2]i and the 
abundance of H is a factor of 10 higher than for He. 

As discussed in Sect. 2.2.1, the threshold energy for carbon 
ejection via nuclear excitation is not well-constrained. We con- 
sider To = 7.5 eV to be a reasonable value, but experimental 
determinations are necessary. Fig. [9] illustrates how the fraction 
of ejected C-atoms changes as a function of the adopted value 
for the threshold energy. The curves show the cumulative effect 
of H, He and C, calculated for To - 4.5, 7.5 and 15 eV in the in- 
ertial and thermal case. Both in the inertial and thermal case, the 
curves corresponding to the various thresholds follow the same 
trend, and for each value considered of To the inertial destruction 
dominates at low velocity and the thermal destruction at high ve- 
locities. As expected the fraction of ejected C-atoms increases 
for decreasing To in the inertial case, while the curves shift to 
the left in the thermal case, implying that the PAHs will start 
to experience significant damage at lower shock velocities. Our 
results also show that, even assuming a high threshold energy, 
PAHs experience a substantial loss of carbon atoms, which is 
complete for velocities above 175 km s _1 is all cases. 

Finally, we investigated how the nuclear destruction process 
depends on the size of the PAH. Fig. [8] shows the fraction of 
ejected carbon atoms from a big PAH with Nc = 200. The de- 
struction of a 200 C-atom PAH follows the same trends with 
shock velocity as for the 50 C-atom case and the curves are al- 
most identical. This is due to the fact that the velocity and tem- 
perature profiles for the 50 and 200 C-atoms molecules are quite 
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similar, and the collision rate and F^ scale linearly with A^c in 
both the inertiai and thermal case (see Eq. [T7ll2"Tl andl2"2"l>. 

4.2. PAH destruction via electronic interaction by ion 
collisions 

Inspection of Fig. Q and [8] reveals that electronic excitation by 
impacting ions plays only a marginal role in the destruction pro- 
cess. For both PAH sizes, carbon is unimportant because of its 
very low abundance. In the inertiai case H does not contribute 
and He contributes marginally at the highest shock velocities, 
while in the thermal case they lead to a substantial atomic loss 
only for A^c = 50 in the highest velocity shock (200 km s ). For 
a 50 C-atom PAH, the low destruction rate due to electronic ex- 
citation reflects the small cross section for this process for these 
low velocity shocks. The inertiai velocities of the PAH lead to 
electronic excitation only being important for the highest shock 
velocities where the impacting ions have a high enough tempera- 
ture to excite the PAHs sufficiently (cf. MJT). The larger number 
of modes available in 200 C-atom PAHs, makes the electronic 
excitation of such PAHs completely negligible over the full ve- 
locity range of the shocks considered here. 

4.3. PAH destruction due to electron collisions 

The fractional carbon atom loss due to collisions with ther- 
mal electrons is also shown in Fig. [7] and [8] For Nc = 50, the 
number of ejected carbon atoms rises sharply above 75 km s _1 , 
leading to total destruction above 100 km s _1 . For Nc = 200, the 
damage is negligible up to 100 km s _1 , increases significantly 
beyond that and leads to complete destruction above 150 km s _1 . 

The energy transferred by impacting electrons rises sharply 
for velocities in excess of 2xl0 3 km s _1 . This results in a dis- 
sociation probability p shaped as a step function: for v > 2xl0 3 
km s _1 p jumps from values close to zero up to 1. This limiting 
velocity applies to a 50 C-atoms PAH; for Nc = 200 the value 
is higher (4xl0 3 km s _1 ), due to the fact that for a bigger PAH 
more energy has to be transferred for dissociation. These veloc- 
ities correspond to electron temperatures of 10 5 K and 3xl0 5 K, 
which are reached for shock velocities of approximately 100 and 
150 km s _1 , respectively. 

4.4. Summary 

A summary of our findings from Fig. [7] and [8] is presented in Fig. 
[TOl This shows the fractional atomic loss, Fj_, due to electron 
and ion collisions, calculated for the two PAH sizes Nc = 50 and 
200. To show how the fractional loss changes as a function of the 
adopted value for Eq, we added the results obtained assuming 
for the electronic dissociation energy the values 3.65 and 5.6 eV, 
lower and higher respectively than our standard value 4.58 eV, 

For ionic collisions, F^ is determined by nuclear interaction: 
inertiai for low Vs and thermal for high vs. For the 50 C-atoms 
PAH, the electronic contribution emerges for the lowest value 
of Eq, 3.65 eV. As already mentioned, the dissociation probabil- 
ity increases for decreasing Eq, so we are not surprised to find 
Fl enhanced by electronic excitation (both inertiai and thermal). 
The electronic contribution desappears for higher dissociation 
energies, as demonstrated by the coincidence between the curves 
for Eq = 4.58 and 5.6 eV. For A^c = 200, all three ionic curves 
are coincident, indicating that electronic excitation by impacting 
ions does not contribute to PAH destruction below 150 km s _1 . 



Above this value, the carbon loss due to electronic interaction is 
covered by the other processes. 

The shift between the ionic curves for the two PAH sizes is 
due to the small differences in the velocity profiles - due to be- 
tatron acceleration -, which imply in the inertiai case a slightly 
higher damage for the bigger PAH. In the thermal case F^ is 
instead independent on the PAH size because the number of 
ejected C-atoms scales linearly with A^c (see Sect. 4.1). The re- 
sulting effect is an almost linear rise for both PAH sizes up to 
150 km s _1 , beyond which the destruction is complete. 

The 'Electron' curves reproduce the behaviour observed in 
Fig. [7] and [8] The comparison with the ionic F^ clearly indi- 
cates that a 50 carbon atoms PAH is already damaged in a non- 
negligible way in low velocity shocks (50 - 75 km s _1 , ion col- 
lisions) and is totally destroyed above 100 km s _1 (electron col- 
lisions). When Eq = 3.65 eV, the fraction of ejected C-atom for 
vs = 75 km s _1 increases from 0.1 to ~0.35 but the minimum 
shock velocity for complete destruction is unchanged (100 km 
s ). For vs > 75 km s _1 , the 5.6 eV curve is almost parallel 
to the 3.65 eV curve and shifted by 25 km s _I toward higher 
shock velocities. For the 200 C-atom PAH, the carbon atom loss 
is dominated by ionic collisions for shocks with velocity below 
100 km s _1 . Above this value, the combined effect of ions and 
electrons leads to a complete destruction. A lower electronic dis- 
sociation energy shifts back by 25 km s _1 the minimum shock 
velocity required for total carbon ejection. When Eq = 5.6 eV, 
destruction starts to be important only above 125 km s _1 , and 
becomes almost complete at 150 km s _1 . 

We adopt the quantity F^ as destruction efficiency to calcu- 
late the PAH lifetime in shocks. 

4.5. Uncertainties discussion 

The main sources of uncertainties which have to be considered 
for this study are related to the adopted shock profiles, to the 
accuracy of the fitting function for the ZBL nuclear stopping 
cross section, and to the choice of an appropriate value for the 
nuclear threshold energy Tq and for the electronic dissociation 
energy Eq. 

The uncertainties related to the adopted shock profiles here 
are principally due to our assumption that we can equivalently 
treat a small, two-dimensional PAH molecule as a small three- 
dimensional grain. In the calculation of the PAH velocity profiles 
through the shock s we use the same formalism as for the grains 
dJones et al.lll996l) . i.e., we assume that the PAH behaves as a 
three-dimensional grain of the same mass. Any uncertainties are 
then due to the inherent differences in the cross section to mass 
ratios for PAHs and grains. As mentioned in Sect. 3, once the 
PAH cross section is averaged over all possible orientations, the 
differences in the PAH and grain cross sections turn out to be 
only of the order of 11%, for a 50 carbon atom PAH, and are 
therefore rather small compared to the other uncertainties that 
we discuss here. 

The accuracy of the ZBL nuclear stopping cross section de- 
pen ds on the ac c uracy of the single analytical function used 
by IZiegler et al.l (Q985) to calculate the interatomic poten- 
tials between atoms. This universal function has been com- 
pared with experimentally determined potentials, with a re- 
sultin g standard deviation betwe en theory and experiment of 
5% (lO'Connor & Biersack1 [l986). An additional test has been 
made comparing the results from the ZBL function with much 
more complex theoretical calculations including more effects. 
In this case as we ll the results agree within few percent (see 
Ziegle ret al.ll 19851 and references therein). 
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Fig. 7. The fractional C-atom loss Fl, due to collisions with H, He and C ions and electrons, as a function of the shock velocity. 
Fl is defined as the total number of ejected carbon atoms divided by the initial number A^c of C-atoms in the PAH molecule. 
The destructive effect of thermal and inertial 'sputtering' induced by collisions with the projectiles are shown. In the nuclear case, 
each lost atom is the result of a single collision with a given projectile, in other words, every collision removes a C-atom from the 
molecule so, the number of destructive collisions equals the number of ejected carbon atoms. For electrons as projectiles and for 
electronic interaction, each collision leads to the ejection of two carbon atoms. The points are calculated for a 50 C-atoms PAH 
assuming the nuclear threshold energy 7q = 7.5 eV and the electronic dissociation energy Eq = 4.58 eV. 
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The main source of uncertainty in the nuclear stopping cal- 
culation is the choice of the threshold energy To. This quantity 
is not well constrained (see Sect. 2.2.1) so we explored a set of 
plausible values. In Fig. [9] we plotted together the values for 
resulting from the total effect of H, He and C, calculated in both 
inertial and thermal case for three different threshold energies 
To = 4.5, 7.5 and 15 eV. The curves corresponding to the highest 
and lowest threshold r™ x and T™" identify a region which can 
be interpreted as the variation in the amount of destruction due to 
the uncertainty in the threshold energy. For the inertial case, this 
uncertainty introduces a variation in the destruction efficiency of 
a factor less than about 2. For the thermal case, the uncertainty 
introduces a shift of the critical shock velocity above which ther- 
mal destruction is dominant from about 100 to 150 km s _1 for 
To ranging from 4.5 to 15 eV. 

An important issue for the electronic and electron stopping 
calculation is the choice of the value for the parameter Eo- We 
adopt the value 4.58 eV, which has been extrapolated for inter- 
stellar conditions from experimental data. Unfortunately the ex- 
trapolation procedure is very model-dependent, so the same set 
of experimental data can lead to significantly different values 
for the interstellar Eo- The problematic fragment binding energy 
is extensively discussed in MJT, the result of a different choice 
for Eo (3.65, 4.58 and 5.6 eV) on PAH processing by shocks is 
shown in Fig. [10] The differences are quite significant, indicating 
the importance of experimental studies on the critical energy Eo 
describing the dissociation probability of highly excited PAHs. 
To summarize, the errors related to the shock profiles and ZBL 
fitting function are quite small, for To we identified a range of 
plausible values, but experimental determinations would be de- 
sirable, while the choice of the parameter Eq is very uncertain 
and urgently requires a better determination. 



5. Discussion 

5.1. PAH lifetime in shocks 

To calculate the timescale for supernova shock waves to destroy 
the interstellar PAHs in the Galaxy, ?snr, we adopt the sam e 
approach used in our previous works ([Jones et al.ll 1994 [l996), 
which is based on the method of McKeeT Tl989l) : 



fSNR 



M ISM 



(l/T SN )/ e (v s )dM s (v s ) 



(28) 



where Mism = 4.5xlO 9 M0 is the mass of the Galactic interstel- 
lar medium (gas and dust including P AHs), t! n = 1 25 yr is the 
effective interval between supernovae (McKee 1989),e(vs) is in 
this case the efficiency of PAH destruction by a shock of velocity 
vs and Ms is the mass of gas shocked to at least vs by a super- 
nova remnant in the Sedov-Taylor stage (in this stage the energy 
is conserved, so that Ms <x E). In a three-p hase model of the 
interstellar medium (Mc Kee & Ostrikei1ll977l) . with a ratio be- 
tween warm and hot intercloud filling factor of / w //h = 0.3/0.7 
= 0.43, the mass Ms and the timescale ?snr become 



M s (vs) = 2914 M Q /v 2 sl 



fSNR - 



9.7 x 10 7 



j e(vs7)/vs 7 dv S 7 



yr 



(29) 



(30) 



were vs7 is shock velocity in units of 100 km s 1 and the assumed 
mean supernova energy is 10 51 erg. 



Using our calculated fractional destruction data, we derived 
analytical expressions for the destruction efficiency e(vs7) for 
electrons and ions. The ionic term represents the total contri- 
bution of all considered ions (H, He and C) and processes (nu- 
clear and electronic stopping, inertial and thermal). Below 150 
km s _1 , we adopt a power law fit for Nc = 50 and a linear fit for 
Nc = 200, while for the remaining velocities the efficiency is 1 . 
For electrons, for F^ < 1, the destruction efficiencies are well 
fit by two linear functions (Nc = 50) and a logarithmic function 
(Nc = 200). The analytical fits reproduce the calculated data 
within few percent. The functional form, fitting parameters and 
corresponding timescales calculated from Eq.[30]are reported in 
Table El 

For electron collisions we find ?snr = 1.6 xlO 8 yr for A^c = 
50 and 4.0 xlO 8 yr for Nc = 200. In case of ion collisions 
the lifetimes are 1.8 xlO 8 yr and 1.4 xlO 8 yr for the small and 
big molecule respectively. The largest uncertainty in these life- 
times result from the uncertainty in the values adopted for To 
and Eq (cf. Sect. 4.5). If we assume the lowest values considered 
(4.5 and 3.65 eV), these lifetimes decrease to 9.2xl0 7 yr and 
8xl0 7 yr (for Nc = 50 and 200 respectively), while the maxi- 
mum values considered for these energies result in lifetimes of 
2.5x10 s yr and 3.3xl0 8 yr. for the small and big PAH. Thus, 
while these parameters are quite uncertain, the derived values for 
the lifetimes are quite robust. Essentially, PAHs are destroyed by 
shocks larger than about 100 km/s and, typically, interstellar gas 
encou nters such shocks once every 100 million years dMcKed 
Il989h . 

From the results obtained assuming our standard values (To 
= 7.5 eV and Eo = 4.6 eV), we argue then that small PAHs are 
preferentially destroyed by electrons whereas big PAHs are more 
affected by ions. 

Our derived values for ?snr for PAH s are signi fi cantly 
shorter than the 6 xlO 8 yr calculated by iJones et al] ([l996) 
for graphite/amorphous carbon grains in the warm intercloud 
medium. Ignoring betatron acceleration, the total number of col- 
lisions per C-atom required for stopping the inertial motion is in- 
dependent of grain size. The difference in lifetimes reflects then 
a difference in sputtering efficiency. This is not surprising be- 
cause of the different approaches adopted for the ion-particle in- 
teractions for PAHs and grains. In both cases, and for the shocks 
that we consider here, it is the sputtering processes that com- 
pletely dominate dust destruction. In the case of grains not every 
atom that is 'knocked on' by an incident ion is lost, i.e. sputtered, 
from the grain. The displaced atom is often embedded deeper 
into the grain and therefore not sputtered from the grain even 
when the displacement energy significantly exceeds the thresh- 
old energy for target atom displacement. This is reflected in the 
fact that, for grains, the sputtering yield is usually much less that 
unity. The shorter lifetime for the PAHs is then be ascribed to the 
fact that in any incident ion interaction the target atom is always 
lost from the PAH when the energy to displace it is greater that 
the required threshold energy. In this case the equivalent PAH 
sputtering yield is then unity. 

Our derived PAH lifetime is much closer to the value of ~2x 
10 8 yr found for a size distribution of hydrogenated amorphous 
ca rbon (a-C:H) grains, typical of th e diffuse interstellar medium, 
by Serra Dfaz-Cano & Jones (2008). This is perhaps just coinci- 
dental; it has its origin in lower average binding energy of C- 
atoms in amorphous carbon than in graphite. a-C:H grains are 
found, as we find for PAHs, to be more susceptible to sputter- 
ing er osion than graphite/amorphous carbon grains (IJones et al .1 
1996) . However PAHs are much more susceptible to erosion than 
a-C:H in fast shocks (> 150 km s~') and this is simply due to 
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Fig. 9. Carbon atom loss due to collisions with thermal and inertial ions (H + He + C) via nuclear interaction. The figure shows the 
comparison betwen the three threshold values To = 4.5, 7.5 and 15 eV for a 50 carbon atom PAH. 

Table 3. PAH destruction analytical fit parameters and survival timescales for electron and ion collisions. 





Electron 


Ion 






Fit 






N c = 50 


N c = 200 


N c = 50 


N c = 200 




parameters 




e(v S 7) = mix + qi 


0.50 < v S7 < 0.75 








mi 


= 0.464 


9i = 


- 0.232 


s(v S7 ) = m 2 x + q 2 


0.75 < v S7 < 1.00 










= 3.536 


12 = 


- 2.536 


e(vsi) = m 3 x + q 3 








0.50 < v S7 < 1.50 


m 


= 0.879 


13 = 


-0.319 


e(v S7 ) = ax b 






0.50 < vs7 < 1.50 




a 


= 0.405 


b = 


2.232 


s(v S7 ) = c \n(x?)/x 




0.50 < v S7 < 1.50 






c 


= 5.366 


d 


= 0.7 


s(v S7 ) = k 


1.00 < v S7 < 2.00 


1.50 < v S7 < 2.00 


1.50 < vs? < 2.00 


1.50 < v S7 < 2.00 




k 


= 1 




«snr (yr) 


1.6 xlO 8 


4.0 xlO 8 


1.8 xlO 8 


1.4 xlO 8 











Note: Fit parameters and survival timescales calculated assuming To = 7.5 eV and Eo = 4.6 eV. 



their small sizes. Thermal sputtering in the hot post-shock gas of 
fast shocks is proportional to the surface area and, small particles 
having a larger surface area per unit mass than large particles, are 
more rapidly eroded. 

5.2. Astrophysical implications 

As shown in Sect. 5.1, the PAH lifetime against shock destruc- 
tion is much shorter than the Stardust injection timescale into the 
interstellar medium f; n j = 2.5 xlO 9 yr. Gas shocked to velocities 
of the order of 50 - 150 km s~' is observed in many regions of 
the in terstellar medium: e.g. toward the star £ Ori (Welt v~et~al1 
|2002|) . in Herbig-Haro jets in the Orion and Vela star forming 
regions ( Podio et"al] 120061) . and in the local interstellar cloud 
(Slavinl l200C Hence, according to our calculations, PAH de- 
struction should be widespread in the ISM. 

If we assume globally that the same f; n j holds for PAHs as for 
'standard' interstellar dust, i.e., that dust and PAH formation are 
coeval in and around evolved stars, our calculated survival times 
for PAHs, i.e., ?snr, indicate that they need to be re-formed in 
the ISM even more rapidly than the larger interstellar grains. If 
PAHs are formed by the fragmentation of larger carbonaceous 



grains then their 'effective' survival time must just be the same as 
that of the larger grains from which they originate. However, this 
can only be true for low velocity shocks, or turbulent regions of 
the ISM, where grain-grain collisions at relatively low velocities 
(of the order of a few km s~') can form PAHs via fragmentation 
and where there is no associated destructive process in operation. 

On a region by region basis our results indicate that PAHs 
should not exist in environments shocked to high velocities (> 
100 km s~'). We conclude that PAHs that exist in unshocked re- 
gions do not survive the passage of shocks with velocities above 
100 - 150 km s~' (depending on their size). They are in fact de- 
stroyed rather early in in the shock at shocked column densities 
of the order of 10 16 - 10 17 cm . Any 'daughter' PAHs produced 
in the post-shock region, by grain fragmentation in grain-grain 
collisions at shocked column densities of the order of 10 17 - 10 18 
cm 2 , will be destroyed by erosion due to their high injection ve- 
locities into the gas following the fragmentation of t heir larger, 
'pare nt' grains that undergo betatron-acceleration (iJones et al.l 
1996). Thus, high velocity shocks destroy all the PAHs that they 
both interact with and produce by fragmentation in high velocity 
grain-grain collisions. 
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Fig. 10. The fractional C-atom loss F^ as a function of the shock velocity, calculated for three values of the parameter Eq: 3.65, 
4.58 and 5.6 eV. F^ is defined as the total number of ejected carbon atoms divided by the initial number Nc of C-atoms in the 
PAH molecule. For the two PAH sizes Nc =50 (top panel) and A^c = 200 (bottom panel), the curves labelled Ton' illustrate the 
cumulative effect of all destructive processes for the ions considered in this study: inertial and thermal 'sputtering' due to nuclear 
and electronic excitation during PAH collisions with H, He and C in the shocks. The curves labelled 'Electron' show the destructive 
effect of thermal electrons. We assume Tq = 7.5 eV. 



In contrast observations show that PAHs lock up about 
3% of the elemental carbon in the ISM (cf. iTielensl |2008l) . 
Disentangling these two scenarios is not easy observationally. 
Here we now consider the case where the PAH emission is as- 
sumed to come from within the shocked region. This scenario 
requires an efficient (re)formation route for PAHs in the dif- 
fuse ISM (see above). However, this is difficult to understand 
since PAHs are a product of high temperature chemistry involv- 
ing abundant carbon bearing precursors such as CH4 and C2H2. 
In the low temperature diffuse ISM, an O-rich environment, 
these precursor species are never really very abundant. This co- 
nundrum is very reminiscent of the general problem of rapid 
dust destruction in the ISM and the long injection time scale 
for freshly formed dust dDwek & Scaloll980tlDraine & Salpeter 
119791: Uones et al.lll994 119961) . We note that, while dust grains 



may be rapidly covered by (t hin)protective coatings between 
successive shock passages (cf. lTielens|[l998l) . this is not a way 
out of this conundrum for PAHs. Perhaps, PAHs can be formed 
through prolonged photolysis of ice mantles accreted inside 
dense molecular clouds, although, presently, there is no exper- 
imental support for this suggestion. Lastly, observations have re- 
vealed the presence of PAHs associated with hot shocked gas 
in stellar (eg., M17) and galactic (e.g., M82) wind regions (cf. 
Sect. 1). As will be discussed more extensively in MJT, these 
PAHs likely trace entrained cold gas which has not been fully 
exposed to the destructive effects of high velocity shocks. 

For vs < 75 km s , PAH are disru pted only by inertial 'sput- 
tering ' due to nuclear interactions. ISerra Dfaz-Cano & Jonesl 
(2008) have studied the erosion of hydrogenated amorphous car- 
bon (a-C:H) arising from ion irradiation in shocks in the warm 
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Fig. 11. The evolution of a 50 carbon atom PAH following the loss of A^c(lost) carbon atoms, as a function of the shock velocity, for 
the two limiting cases: instantaneous and random removal of the lost carbon atoms (top row), and carbon atom removal only from 
the periphery of the molecule (bottom row). 



intercloud medium, using exactly the same shock profiles that 
we use for our PAH study, and it is interesting to note that, for 
a 50 carbon atom PAH molecule, the percentage of destruction 
is the same as for a-C:H, indicating that the lower density, more 
easily sputtered a-C:H and PAHs exhibit similar erosion charac- 
teristics in lower velocity shocks (^ 125 km s~'), as discussed in 
the previous section. 

We now consider what happens to a PAH as a result of the 
ejection of aromatic carbon atoms, by the impacting ions, as a 
function of the fraction of carbon atoms removed from the struc- 
ture. Our results show that ionic collisions severely modify PAH 
in shocks with velocities between 75 - 150 km s . Nuclear and 
electronic interaction lead to distinctly different molecular de- 
struction routes. Specifically, electronic excitation (either by im- 
pacting ions or electrons) lead to a high vibrational excitation of 
the PAH and this PAH will relax by 'losing' the weakest link in 
its skeleton. Initially, this will be the peripheral H-atoms or other 
functional groups. For large PAHs (50-100 C-atoms), the result- 
ing 'pure-C may quickly isomerize to very stable carbon cluster 
such as fullerenes. Fullerenes are very stable against C2 loss. 
Measured dissociation energies of fullerenes are i n the range of 
7 to 9 .5 eV, with C 60 itself at E = 9.8 ± 0.1 eV dTomita et all 
1200 lh . 

In contrast, nuclear interaction will act on the C-atom hit by 
the projectile ion. This will randomly remove C-atoms from the 
C-skeleton of the PAH. Some isomerization may occur if inter- 
nal excitation energy is left behind. However, unless the C-atom 
loss in the shock is very large, likely this is insufficient to affect 
the overall PAH structure. In a simple 'geometrical' analysis, if 
one randomly removes C atoms from a PAH with Nq = 50, in the 
absence of any annealing of the PAH chemical structure, we find 
that the loss of 10% of the C atoms leads to a loss of of the order 
of « 50 % of the aromatic character. However, the loss of a; 20 % 
of the C atoms leads to the almost complete loss of aromatic 
character and to the onset of the fragmentation of the molecule. 



Fig.fTTIillustrates these effects and shows the PAH evolution fol- 
lowing the loss of carbon atoms, A^c(lost), for the two limiting 
cases: 1) where there is an instantaneous and random removal of 
the lost carbon atoms (appropriate for nuclear interaction) and 
2) where the carbon atoms are removed only from the periphery 
in order to preserve aromatic domain as much as possible (likely 
appropriate for electronic excitation through either impact ions 
or electrons). The reality of PAH erosion in shocks may well lie 
somewhere between these two extremes and may also involve 
isomerisation and the formation of five-fold carbon rings that 
distort the structure from a perfectly two-dimensional form. This 
then begs the question as to the exact form and structure of small 
carbon species once growth resumes by atom insertion and ad- 
dition. 

After the shock wave has passed, the resulting PAH can re- 
act chemically with impacting H, C, O, and N atoms. A dis- 
tinction can be made between H and O in the one hand and C 
and N on the other hand, The latter can restore the stable aro- 
matic character of the PAH skeleton while the former lead only 
to very unstable structures. Likely, then, H and O atom addi- 
tion can be reversed by U V photon absorption without loss of C 
(lAllamandola et alll989l) . Hence, the PAH carbon skeleton may 
be able to 'repair' itself to an aromatic structure, possibly incor- 
porating N-atoms. The possible presence of N-atoms deeply in 
the C-skeleton of PAHs in the diffuse medium (but not in C-dust 
birth sites such as C-rich planetary nebula e) has been infer red 
from the peak position of the 6.2 /im band dPeeters et al.l2002T) . It 
has also been suggested that the 6.2 //m band position variations 
can be explained by a varying aliphatic to arom atic carbon con- 
tent in carbonaceous particles dPino et a l. 2008). The difference 
in molecular structure between circumstellar PAHs and the gen- 
eral ISM implies an active chemistry in the ISM that is able to in- 
sert N atoms deeply in the carbon skeleton. This is very puzzling 
given the very stable character of the aromatic backbone and the 
low temperature of the ISM. Here, we surmise that, unlike UV 
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photolysis or chemical attack, nuclear interactions in interstellar 
shocks may be a viable pathway to promote N-incorporations 
deep inside PAHs. However, given that the nitrogen abundance 
is a factor of a few lower than that for carbon - depending on the 
fraction of carbon remaining in dust and PAHs - carbon atom in- 
sertion ought to be favoured over that for nitrogen, in the absence 
of any chemically-selective route for nitrogen insertion. 



We surmise that the molecular structure of PAHs is 
strongly affected by shock processing in the interstellar medium. 
Electronic excitation by impacting ions or electrons may lead 
to isomerization into stable pure-C species such as fullerenes. 
In contrast, nuclear interaction may lead to the formation of 
N-containing PAHs. Further laboratory studies are required to 
demonstrate the viability of these chemical routes. 



6. Conclusions 

We have extensively studied the effects of PAH processing by 
shocks with velocities between 50 and 200 km s _1 , in terms 
of collisions with ions and electrons which can lead to carbon 
atom loss, with a consequent disruption and destruction of the 
molecule. 

An ionic collision consists of two simultaneous processes 
which can be treated separately: a binary collisions between the 
projectile and one single atom in the target (nuclear interaction) 
and the energy loss to the atomic electrons (electronic interac- 
tion). For the nuclear interaction, we modified the existing the- 
ory in order to treat collisions able to transfer energy above a 
specific threshold. This is the case we are interested in, which 
has not been treated in previous studies. For electronic interac- 
tion and collisions with electrons we developed specific models 
for PAH, described in MJT. 

The PAH dynamics in the shoc ks is evaluated using the same 
approach as in our previous work (I Jones et aT I fl994l[l996h . For 
nuclear interaction, the level of carbon atom loss increases for 
decreasing values of the threshold energy Tq. We adopt Tq - 7.5 
eV as a reasonable value, but experimental determinations are 
necessary. In ionic collisions, the carbon contribution to PAH de- 
struction is totally negligible because of its very low abundance 
with respect to H and He. The fractional destruction induced by 
nuclear excitation increases with the PAH size, while in case of 
electronic excitation and electron collisions is lower for higher 
Nq values, i.e. bigger PAHs are more resistant than the smaller 
ones against electron and electronic processing. 

The parameter Eq required for the evaluation of PAH de- 
struction due to electron and electronic interaction is unfortu- 
nately not well constrained. We adopt a value of 4.58 eV con- 
sistent with extrapolations to interstellar conditions, but better 
determinations would be desirable. 

Electronic interaction, both inertial and thermal, plays a 
marginal role in PAH processing by shocks. We find that 50 car- 
bon atoms PAHs are significantly disrupted in ionic collisions 
for shock velocities below 75 km s _1 , mainly by inertial 'sput- 
tering' by helium during nuclear interaction. Our results indicate 
5- 15 % C atom loss, sufficient to cause a severe de-naturation of 
the PAH aromatic structure. Above 100 km s _1 such PAHs are 
instead totally destroyed by collisions with thermal electrons. 
For Nc = 200, PAHs experience increasing damaging caused 
by nuclear ionic collisions up to 100 km s _1 , which turns into 
complete atomic loss for higher velocities. In this case the de- 
struction is due to the combined effect of electrons and nuclear 
interaction with thermal ions. 

The calculated PAH lifetime against destruction, fsNR> is 1 -6 
Xl0 8 yr and 1.4 xlO 8 yr for N c = 50 and 200 respectively. Small 
PAHs are preferentially destroyed by electrons, big PAHs by 
ions. The calculated lifetimes are smaller than the values found 
for carbonaceous grains (6 xlO 8 yr) but close to that for hydro- 
genated amorphous carbon ( 2 XlO 8 yr), and far from the star- 
dust injection timescale of 2.5 xlO 9 yr. The presence of PAHs 
in shocked regions therefore requires an efficient reformation 
mechanism and/or a protective environment. 



Appendix A: S n , cr and (T) 

A fundamental quantity to describe the nuclear scattering is the 
energy transfer cross section cr(E, T), which is function of the 
kinetic energy E of the projectile and of the energy T transferred 
to the target by the projectile in a single collision. At low en- 
ergies (e < 1), an approximated expression for the cross section 
can be calculated using the power approximation to the Thomas- 
Fermi model of interatomic interaction, i.e. with a potential of 
the form V(r) cc r~ l l"\ where r is the distance between collid- 
ing nuclei and m is a parameter related to the steepness of the 
interatomic potential. The quantity m can also be interpreted as 
an indicator of the energy of the projectile, varyi ng slowly from 
m = 1 at high energies to m » at low energies ( Lindhar d et alj 
1968; Winterb on et ~il l970). We have that 



do-(E, T) =* C,„ E~ m r 1 -'" dT < T < T, 



with 



n 



Mi 



2Z x Zie 2 



2m 



(A.l) 



(A.2) 



where T m is the maximum transferable energy, corresponding to 
a head-on collision (impact parameter p - 0). The dimensionless 
quantity A m = A m (m) varies slowly fro m 0.5 for m _ 1 (high 
energy, i.e. pure Rutherford scattering, I Simmon IClll, to 24 
for m — (very low energy). 

The power approximation of the Thomas-Fermi cross sec- 
tion (Eq. I A. 1 b gives the following expression for the nuclear 
stopping cross section, obtained evaluating the integral in Eq. |4] 
between and T m : 



S n (E) = 



1 



1 — m 



1-m 77l-2m 



(A.3) 



For Thomas-Fermi interaction, using Eq.rT HA~2l l3land rA.3l this 
leads to 



S n (E) = 47raZ l Z 2 e 2 
with 

■Sn(e) 



Mi 



Mi + Mi 



SnO) 



, 1 -2m 



2(1 -m) 



(A.4) 



(A.5) 



For heavy screening (s«l,ras 0) the accuracy of the Thomas- 
Fermi reduced stopping cross section s n {s) to reproduce the ex- 
perimental data is at best a factor of two. For this reason we adopt 
inste ad the Univer s al red uced stopping cross section s^ (Eq. [5]) 
from lZiegler et alJ d 19851) with the appropriate screening length 
«u (Eq. [2j, which provides a good fit to the experimental data at 
low energies as well. 

The nuclear stopping cross section above threshold is given 
by Eq. [10] Remembering that T m = yE and Tq = y£o n , this can 
be rewritten as follows 
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Using a similar approach we can write the term (T m m /m) as 



C F 



Sn(E) = -^—[TZr-T'-"] 



r m '» T m m l-m T n 
m m 1 — m T„ 



r F~ m C F~ m 

L ^1-m _ '-m L 

l-m m ~ l-m 



l-m 1 
X X 



1 — m 



1 — m 



c £ 0n 



1 — m m r„ 

\-m E \-m l _ m J 
X X . 

l-m m 7.E 



The first term in the right side of the equation is equal to the 
nuclear stopping cross section in the no-threshold case (To = 0) 



Then second term in the right side of the equation then becomes 



Cm 7 



l-m 



1 — m 



E 1 - 2 '" = AnaZ 1 Z 1 e 2 — ^— s n (e) ^ Sj|(£). 
Mi + M 2 



y-t p-m in 

m l-m 



£ l-2m x 



1 — m 1 



x 



m 7.E 



l-m 1 



^(E)x x 



In the second term we use the equality E "' = (.E 1 2m /£ 1 m ) to 
obtain 

L-m 7 rr-m pl-m _ 7 E>l-2m / '-'On 



m 7Z1 



Combining the two terms we obtain for cr(E) the final expression 



1 — m 



■ E~ m En 



1 — m 



cr(E) = AnaZ\Z 2 e 



2 My 



-s n (s) 



1 — m 1 



Mi + M 2 m yE 



m 



i 



= A 7r a Zi Z 2 e — — s„(e) x — 

Mi + M 2 \ E I 



The expression for the average transferred energy (T(E)} de- 
rives directly from Eq. Q~2] when applying the relations T m = yE 
and T = yE 0n 



Combining the two we obtain the following expression for S n {T(E)) - 



S n (E) = AnaZiZ 2 e 



2 Mi 



Mi + M 2 



-(f)' 



The total cross section is given by Eq. QT| 

C F~ m 



cr(E) 



1^0 J 



m l m - 1 



1 - m T~ m - r m m 



m r 



l-m 



(Pl-m^l-m) 



1 - m y-m ( E -m _ E - m \ 



£?l-m _ j?l—m 



l-m £ ™ - E~ m 



m m 



Cp-m j i e»— m ' m 

m ^ 



Appendix B: Low and high energy regime above 
threshold 



The term (T m /m) can be rewritten as follows 

f-m 71 l-m i _ m f-m 
X X 



m l-m m r,l 



l-mpl-m ,_ m y -« £ -« 
X X 



"On 



1 — m m yl m £\-m 

y \-m E l-m 1 _ m , , E{) 



X X 



1 — m m 7 



yE \ E ) 



Then we have 



In the low energy regime m can be taken equal to 0, so the dif- 
ferential cross section dcr becomes 



dcr(E,T) = s -Aoa 2 T- l dT. 



(B.l) 



The corresponding expressions for S n (E), cr(E) and (T(E)) 
above threshold are calculated from Eq. [7j Eq. [8] and Eq. [9] re- 
spectively 



S a (E) = -Aoa 2 y(E-E 0a ) 



n , E 
cr (E) = — Ana In 



Eon 



(B.2) 
(B.3) 



C F 



Tn'" C m y x m „ 1 _ 2m ., 1 -m _ 1 /£()„ 1 
m l-m m y i 



-m 



. £ l-2m x 



X 



1 /£0n\~ 

7£ I £ / 



nV 7 m yE \ E I 



(T(E)) = 7(£-£on) | In — 

^0n 



In the high enery regime m = 1 , then we have 



M 



dcr(E, T)^—Z} Zl e A n- T 1 dT 
M 2 1 2 £ 



1 



(B.4) 



(B.5) 
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and consequently 



,™ Mi 2 2 4 1 1 

cr (£) = Zf Zte n 

M 2 1 2 £ r £on 



<r(£)> = r £ 0n In — 



(B.6) 

(B.7) 

(B.8) 
(B.9) 



Appendix C: Orientation correction 

To calculate the orientation correction factor (S ), let us consider 
the two versors u and n, perpendicular to the PAH surface and 
forming between each other the angle §. This configuration de- 
fines two possible orientations for the molecule. The average ori- 
entation is then given by the following integral: 



(S) 



1 f 
= — un 

2n J 



(CI) 



where dQ = sin#d#d^ is the solid angle element and 
u ■ n - cos & is the scalar product between the two versors. 
With these substitutions we obtain 



(S) 



i r 27r r n< 

Z7r J(p=0 Ji?=i 



OT /2 1 

cos # sin & di? = - 

o 2 



(C2) 
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